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ABSTRACT 

We analyse the Fundamental Plane (FP) of 26 strong gravitational lens galaxies 
with redshifts up to z ~ 1, using tabulated data from Kochanek et al. (2000) and 
Rusin et al. (2003). The lens galaxies effectively form a mass-selected sample of early- 
type galaxies in environments of relatively low density. We analyse the FP and its 
evolution in the restframe Johnson B and Gunn r bands. Assuming that early-type 
galaxies are a homologous family, the FP then provides a direct measurement of the 
M/L ratio evolution. 

If we assume that the M/L ratios of field early-type galaxies evolve as power- 
laws, we find for the lens galaxies an evolution rate dlog(M/L)/dz of —0.62 ± 0.13 in 
restframe B and —0.47 ±0.11 in restframe Gunn r for a flat cosmology with LIm = 0.3 
and 12 a = 0.7. For a Salpeter (1955) Initial Mass Function and Solar metallicity these 
results correspond to mean stellar formation redshifts of (z*) = 1.8+(A5 an d 
respectively. After correction for maximum progenitor bias, van Dokkurn & Franx 
(2001) find a mean stellar formation redshift for cluster galaxies of (zf) = 2 .(Fq' 2, 
which is not significantly different from that found for the lens galaxies. However, if 
we impose the constraint that lens and cluster galaxies that are of the same age have 
equal M/L ratios and we do not correct for progenitor bias, the difference is significant 
and we find that the stellar populations of the lens galaxies are 10-15 % younger than 
those of the cluster galaxies. 

We find that both the M/L ratios as well as the restframe colors of the lens galaxies 
show significant scatter. About half of the lens galaxies are consistent with an old 
cluster-like stellar populations, but the other galaxies are bluer and best fit by single 
burst models with younger stellar formation redshifts as low as z* ~ 1. Moreover, 
the scatter in color is correlated with the scatter in M/L ratio. We interpret this 
as evidence of a significant age spread among the stellar populations of lens galaxies, 
whereas the ages of the stellar populations of the cluster galaxies are well approximated 
by a single formation epoch. 

Key words: galaxies: elliptical and lenticular, cD galaxies: evolution - galaxies: 
formation - gravitational lensing 


1 INTRODUCTION 

One of the central problems in astronomy is that of galaxy 
formation and evolution: when were the visible parts of 
galaxies assembled and when were the stars formed? The 
merging of galaxies leads to changes in their masses, and 
stellar evolution changes their luminosities. The evolution 
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of the mass-to-light (M/L) ratio relates the mass and lumi¬ 
nosity evolution. 

Galaxy mass measurements are notoriously difficult. 
Fortunately, empirical relations such as the Tully-Fisher re¬ 
lation for spiral galaxies (Tully & Fisher 1977) and the Fun¬ 
damental Plane (FP) for early-type galaxies (Dressier et al. 
1987; Djorgovski & Davis 1987) can provide us with infor¬ 
mation on the masses and mass evolution of galaxies. The 
FP is a tight relation between the structural parameters and 
velocity dispersion, which under the assumption of homol- 
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ogy implies that M/L oc M * with low scatter (Faber et al. 
1987). 

Due to stellar evolution the M/L ratio of a stellar pop¬ 
ulation changes with redshift, and hence also the FP will 
change. The redshift evolution of the intercept of the FP is 
proportional to the evolution of the mean M/L ratio. Hence, 
the tightness of the FP relation makes it a very sensitive 
indicator of the mean age of the stellar population of early- 
type galaxies (van Dokkum & Franx 1996). 

For cluster galaxies, the M/L ratio evolves very slowly, 
indicating that the stars were formed at redshifts z > 3 (e.g. 
Kelson et al. 1997; Bender et al. 1998; van Dokkum et al. 
1998). Current semi-analytical hierarchical models place the 
assembly time of typical early-type galaxies at much lower 
redshifts (e.g. Kauffmann 1996; Kauffmann & Chariot 1998; 
Diaferio et al. 2001). However, in hierarchical galaxy forma¬ 
tion models the assembly time and star formation epoch are 
strongly dependent on the environment, with cluster early- 
type galaxies forming much earlier than those in the general 
field (e.g. Kauffmann 1996). Hence cluster early-type galax¬ 
ies do not provide the most stringent tests of these models. 
Moreover, the measured evolution of early-type galaxies may 
underestimate the true evolution because of the effects of 
morphological evolution. If many early-type galaxies evolved 
from late-type galaxies at z < 1, the sample of early-type 
galaxies at high redshift is only a subset of all progenitors 
of present-day early-type galaxies. This would lead us to 
underestimate the luminosity evolution, and hence overes¬ 
timate the stellar formation redshift (see van Dokkum & 
Franx 2001, hereafter vDFOl). 

Recent studies have started to explore the FP and M/L 
ratios of high redshift early-type galaxies in the general field. 
Interestingly, the FP of field early-type galaxies appears to 
be quite similar to that of cluster early-type galaxies out 
to z 0.55 (Treu et al. 2001; van Dokkum et al. 2001), in 
apparent conflict with current semi-analytical models (see 
van Dokkum et al. 2001). However, there are indications for 
a significant offset between the two populations at higher 
redshift (Treu et al. 2002). 

Kochanek et al. (2000, hereafter K00) analysed the FP 
and color evolution of strong gravitational lens galaxies up 
to z ~ 1. The lensing cross section is dominated by galaxies 
with high central velocity dispersions, and the lens galax¬ 
ies effectively provide a mass-selected sample of field early- 
type galaxies which can be compared to optically selected 
samples of field and cluster galaxies. The mass-selection is 
important because it limits Malmquist-type biases, and it is 
much less sensitive to selection effects caused by morpholog¬ 
ical evolution. K00 find that the stars constituting the lens 
galaxies must have formed at z > 2 for a flat cosmology with 
Dm = 0.3 and Ha = 0.7, and conclude from their analysis 
that there are no significant differences between field and 
cluster early-type galaxies. 

K00 analyse the evolution of the FP in the observed 
photometric bands, and their modeling necessarily includes 
the large variation in observed magnitudes and colors due 
to redshift (the “K”-correction). This complicates the mea¬ 
surement of the smaller effect due to evolution of (the stellar 
populations of) the lens galaxies. Rusin et al. (2003, here¬ 
after R03) study the same sample of lens galaxies, with re¬ 
cent photometric observations included. They convert the 
data from observed filters into magnitudes in standard rest- 


frame bands, so that the evolution of the M/L ratio that 
follows from the FP (under the usual assumption that early- 
type galaxies are a homologous family) can be investigated 
instead of the evolution within the FP. This clarifies the 
analysis and allows for direct comparison of the results with 
those of other FP studies. From their M/L evolution anal¬ 
ysis R03 find a (2a) lower limit z > 1.8 for the formation of 
the stars in lens galaxies. Although this limit is more pre¬ 
cise, it is similar to that of K00. Hence, R03 also conclude 
that the evolution measurements favor old stellar popula¬ 
tions among field galaxies, like those of cluster galaxies, and 
argue against significant episodes of star formation at z < 1, 
as predicted by the semi-analytical hierarchical models. 

For the analysis in this paper we use the tabulated data 
of K00, extended with the recent photometric measurements 
as given by R03. We convert the observed surface bright¬ 
nesses and colors to the restframe Johnson B and Gunn r 
bands. Our transformation to restframe differs from that of 
R03. In both approaches the modeled color between rest¬ 
frame band and observed filter is used to convert the ob¬ 
served magnitude into an estimate of the restframe mag¬ 
nitude. R03 obtain the modeled color for a given spectral 
energy distribution, whereas we use the observed color be¬ 
tween a pair of filters to interpolate between the modeled 
colors for four different spectral types (Fig.©. Moreover, 
while R03 use the (weighted) contribution of all observed 
magnitudes to determine the restframe magnitude, we se¬ 
lect the best one with small observational error and filter 
close to the (redshifted) restframe band. Hence, the model 
dependence of our correction is small and we exclude ob¬ 
served magnitudes with large uncertainties, minimizing the 
error in the resulting restframe magnitude. We analyze the 
M/L evolution of the lens galaxies and compare our results 
with those of R03 and with results from previous studies. 
Additionally, we test whether age differences between the 
lens galaxies are significant. To this end, we study the scat¬ 
ter in both the M/L evolution and the restframe colors of 
the lens galaxies, and investigate whether the deviations are 
correlated. 

This paper is organized as follows. In Section 12] we de¬ 
termine the FP parameters, using the lensing geometry to 
estimate the velocity dispersion. The transformation from 
observed to restframe bands is described in Section© In Sec¬ 
tion 2] we construct the FP of the lens galaxies and present 
the M/L evolution derived from the FP. We estimate the age 
of the stellar populations of the lens galaxies by fitting single 
burst models to the M/L evolution in Section © In Section 
©we study the colors of the lens galaxies. The results are 
summarized and discussed in Section© Unless stated other¬ 
wise we assume Hq = 50kms -1 Mpc -1 and a flat cosmology 
with 9 m = 0.3 and Ha = 0.7. We note that our results are 
not dependent on the value of the Hubble constant. 


2 FP PARAMETERS 

The study of the FP of strong gravitational lens galaxies 
differs in two important aspects from that of cluster galaxies. 
First, the lens galaxies are individual galaxies spread over a 
large range in redshifts, instead of an ensemble of galaxies at 
the same redshift. Additionally, whereas studies of the FP 
of optically selected galaxies measure velocity dispersions 
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Name Lens 

z 

cr c 
(km s" 

- 1 ) 

&C* 

(km s _1 ) 

(kpc) 

Ve,B z 

(mag arcsec -2 ) 

M e,r z 

(mag arcsec -2 ) 

0047-2808 

0.49 

254 

± 

26 

229 ± 15 

7.7 ± 0.7 

21.79 

± 

0.16 

21.11 

± 

0.28 

Q0142-100 

0.49 

224 

± 

22 


4.3 ± 0.2 

20.42 

± 

0.05 

19.64 

± 

0.05 

MG0414+0534 

0.96 

303 

± 

30 


8.6 ± 1.6 

21.07 

± 

0.17 

19.97 

± 

0.16 

B0712+472 

0.41 

181 

± 

18 


2.8 ± 0.4 

20.67 

± 

0.17 

19.60 

± 

0.16 

RXJ0911-I-0551 

0.77 

260 

± 

26 


7.0 ± 0.6 

21.74 

± 

0.17 

20.75 

± 

0.14 

FBQ0951+2635 

0.24 

128 

± 

13 


0.9 ± 0.2 

19.82 

± 

0.25 

19.05 

± 

0.25 

BRI0952-0115 

0.41 

117 

± 

12 


0.8 ± 0.2 

20.00 

± 

0.21 

18.61 

± 

0.20 

Q0957+561 

0.36 

431 

± 

43 

305 ± 11 

14.1 ± 1.3 

22.48 

± 

0.12 

21.45 

± 

0.11 

LBQS1009-0252 

0.88 

198 

± 

21 


1.9 ± 0.3 

19.81 

± 

0.14 

18.90 

± 

0.13 

Q1017-207 

0.78 

151 

± 

16 


3.1 ± 0.1 

21.28 

± 

0.49 

20.28 

± 

0.07 

FSC10214+4724 

0.75 

241 

± 

26 


11.8 ± 5.2 

22.97 

± 

0.44 

22.03 

± 

0.56 

B1030+074 

0.60 

218 

± 

22 


4.2 ± 0.6 

21.38 

± 

0.14 

20.14 

± 

0.25 

HE1104-1805 

0.73 

316 

± 

32 


6.4 ± 1.9 

21.34 

± 

0.33 

20.22 

± 

0.30 

PG1115+080 

0.31 

210 

± 

21 

288 ± 27 

3.0 ± 0.1 

21.23 

± 

0.06 

20.01 

± 

0.05 

HST14113+5211 

0.46 

190 

± 

19 


3.8 ± 0.4 

21.70 

± 

0.09 

20.72 

± 

0.11 

HST14176+5226 

0.81 

292 

± 

29 

230 ± 14 

7.5 ± 0.9 

20.98 

± 

0.15 

20.35 

± 

0.12 

B1422+231 

0.34 

160 

± 

16 


2.1 ± 0.6 

21.18 

± 

0.25 

19.93 

± 

0.24 

SBS1520+530 

0.72 

220 

± 

22 


3.5 ± 0.3 

20.13 

± 

0.17 

19.36 

± 

0.08 

MG1549+3047 

0.11 

188 

± 

19 

242 ± 20 

2.3 ± 0.2 

21.39 

± 

0.09 

20.12 

± 

0.09 

B1608+656 

0.63 

292 

± 

29 


6.2 ± 1.0 

20.68 

± 

0.19 

19.88 

± 

0.15 

MG1654+1346 

0.25 

206 

± 

21 


4.9 ± 0.1 

21.99 

± 

0.07 

20.73 

± 

0.05 

MG2016+112 

1.00 

299 

± 

30 

328 ± 32 

2.5 ± 0.3 

19.22 

± 

0.10 

18.06 

± 

0.10 

B2045+265 

0.87 

378 

± 

38 


4.1 ± 1.3 

20.50 

± 

0.44 

19.41 

± 

0.30 

HE2149-2745 

0.50 

203 

± 

20 


4.3 ± 0.4 

20.92 

± 

0.12 

20.36 

± 

0.12 

Q2237+030 

0.04 

168 

± 

17 

220 ± 31 

4.3 ± 0.8 

22.37 

± 

0.50 

21.15 

± 

0.23 

HS0818+1227 

0.39 

251 

± 

25 


6.6 ± 0.2 

22.18 

± 

0.07 

21.13 

± 

0.04 


Table 1. FP parameters of 26 strong gravitational lens galaxies with redshifts up to z ~ 1. For 5 lens galaxies the redshift is not known 
spectroscopically, and a photometrically estimated value (in italics) is given. The velocity dispersion cr c , within the standard aperture 
with a diameter of 3”4 at the distance of the Coma cluster, follows from the lensing geometry, assuming a singular isothermal sphere 
mass model and a Hernquist (1990) luminosity profile. The velocity dispersion from stellar kinematics a * has been measured for 7 lens 
galaxies (references in text). The effective radius r e and effective surface brightness p e follow from fits to an r 1 / 4 law. To allow a direct 
comparison with the local FP, the effective surface brightness has been corrected to the restframe Johnson B and Gunn r band by 
interpolating between filters. 


from spectra, for lens galaxies we use the lensing geometry 
to estimate this quantity. The separation between lensed 
images of background sources increases with the mass of the 
lens and is therefore a measure of the velocity dispersion. 
The two remaining FP parameters, the effective radius and 
surface brightness, are determined from surface photometry 
as for cluster galaxies. 

2.1 Velocity dispersion 

For a singular isothermal sphere (SIS) mass model the re¬ 
lation between the velocity dispersion of the matter dis¬ 
tribution and the separation of the source images Ad is 
Ad = d/ c) 2 Dns/Dos. Here, Dls and Dos are the an¬ 

gular diameter distances from the lens galaxy to the source 
and from the observer to the source. To determine these 
values, the redshifts of both the lens galaxy and the source 
are needed (e.g. Hogg 2000). If the redshift of a lens galaxy 
is not known spectroscopically, we adopt a photometrically 
estimated value 1 with a 10% uncertainty. The lens systems 
for which no redshift is known for the source, we exclude 
from our analysis. This leaves a total of 26 lens galaxies, of 

1 Obtained from the CfA-Arizona Space Telescope Lens Survey 
(CASTLES) web site at http://cfa-www.harvard.edu/castles/ 


which 5 have a redshift that is estimated photometrically 
(Table0. The velocity dispersion an depends also (weakly) 
on cosmology through the angular diameter distances, but 
is independent of the value of the Hubble constant since the 
distances appear as a ratio. 

The velocity dispersion an is that of the total matter 
distribution, including possible dark matter, rather than the 
central velocity dispersion of the stellar component a c . For a 
given mass and luminosity distribution, the ratio g = a c /an 
can be modeled by solving the single Jeans equation for 
a spherical system with a$ = ag and /3 = 1 — a 2 /a 2 the 
anisotropy parameter (e.g. Binney & Tremaine 1987). The 
overall mass distribution is assumed to be isothermal. We 
model the luminosity distribution by a Hernquist (1990) 
profile with characteristic radius a = r e /1.8153. For each 
lens galaxy with effective radius r e in kpc (see Section E3>. 
we integrate the resulting line of sight velocity dispersion 
within the Coma aperture to obtain an estimate for g, so 
that the stellar velocity dispersion follows as a c = g an- 

For the sample of 26 lens galaxies, an is distributed with 
a (biweight 2 ) mean of 230kins -1 and a (biweight) standard 

2 Throughout this paper we use the biweight location and scale as 
estimators of the mean and standard deviation (rms) respectively. 
These estimators are robust for a broad range of non-Gaussian 
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deviation of 61kms _1 . In the case of an isotropic system 
(/3 = 0), the estimated stellar velocity dispersion a c has a 
mean of 223 kms -1 and a standard deviation of 70kms -1 . 
Hence, the properties of the sample of lens galaxies will be 
typical for early-type field galaxies close to L*, which have 
a characteristic velocity dispersion around 225kms -1 (e.g. 
Kochanek 1994, 1996). 

Several sources contribute to the uncertainty in a c . We 
account for the scatter in the image separation Ad, which is 
small and set to 2%. Also the error in the effective radius r e 
(Table0 contributes through the modeled ratio g. For those 
galaxies for which only a photometric estimate of the lens 
redshift zi is available, we included an error of 10 %. Since 
the angular diameter distance depends on redshift this error 
is an additional contribution to the uncertainty in a c via the 
determination of ctd and via the conversion of the effective 
radius in arcseconds into physical units of kpc. Apart from 
these errors in observational parameters, we also want to 
take into account that we have made several assumptions 
in the modeling of the ratio g. For the 26 lens galaxies, the 
distribution of g has a mean of 0.90, 0.94 and 1.01 for an 
anisotropy parameter (3 of —0.5, 0 and 0.5 respectively. For 
each lens galaxy we apply the isotropic case (3 = 0, but to 
take the variation of g with (3 into account, together with 
the assumption of a SIS mass model and a Hernquist profile, 
we assume an additional error of 10 % in g and hence in a c . 

For 7 out of the 26 lens galaxies presented in this paper, 
a velocity dispersion measured from stellar kinematics cr c + is 
available: 0047-2808 (Koopmans & Treu 2003), Q0957+561 
(Tonry & Franx 1999; Falco et al. 1997; Rhee 1991), 
PG1115+080 (Tonry 1998), HST14176+5226 (Ohyama et 
al. 2002), MG1549+3047 (Lehar et al. 1996), MG2016+112 
(Koopmans & Treu 2002) and Q2237+030 (Foltz et al. 
1992). Comparing the measured velocity dispersions (Table 
0 with the modeled a c from the lensing geometry, we find 
for these 7 galaxies that the ratio <j c +/&c is distributed with a 
mean of 1.07 and dispersion of 0.27. For the sample of 7 lens 
galaxies both methods are consistent. However, due to pe¬ 
culiarities of the lensing system, such as the contamination 
of Q0957+561 by the mass distribution of the underlying 
cluster, the velocity dispersion from both methods can be 
significantly different for individual lens galaxies. 


2.2 Effective radius and surface brightness 

The lensing systems were observed with the WFPC2, NICM0S1 
and NICM0S2 camera on the HST, in filters ranging from the 
visual F555W through the infrared F205W filter. For each 
system, the image with optimal contrast between the lens 
galaxy and the images of the background source is selected, 
and the effective radius r e and the mean surface brightness 
within the effective radius (/r e ) from fits to an r 1//4 law is 
determined. The data and model fits are described in detail 
by Lehar et al. 2000 and K00. Note that r e depends slightly 
on cosmology due to the conversion of the effective radius 
from units of arcseconds from the fit into physical units of 
kpc. Moreover, since the lens redshift is used in the con¬ 
version, we have to take into account a (small) additional 


contribution to the uncertainty in r e due to the assumed 
10 % error in the lens redshift in the case it is estimated 
photometrically (see also Section f2.lt . 

In this paper we use the surface brightness at the ef¬ 
fective radius (in magarcsec -2 ), which is related to {fj, e ) 
by fj, e — (g> e ) = 1.393. We also define 7 e = io~ Me / 2 - 5 . In the 
following, we refer to the filter in which the fit was made 
as the reference filter. The dependence of effective radius on 
passband due to color gradients can be ignored because of 
the strong correlation between g, e and r e (see Section 0, 
and we can use the observed colors (tabulated by K00) to 
calculate the effective surface brightness in each filter. 

The effective surface brightnesses and colors are cor¬ 
rected for Galactic extinction with an Rv = A(V)/E(B — 
V) ss 3.1 extinction curve for a diffuse stellar medium (e.g. 
Cardelli, Clayton & Mathis 1989; O’Donnell 1994). The 
galactic extinction E(B — V) is obtained from Schlegel, 
Finkbeiner & Davis (1998). Extinction corrections on the 
effective surface brightness are typically ~ 0.03 magnitudes 
in the reference filter, which in most cases is the F160W 
filter. The three galaxies MG0414+0534, MG2016+112 and 
B2045+265 are exceptions with significantly higher galactic 
extinction of 0.18, 0.14 and 0.14 magnitudes respectively. 


3 TRANSFORMATION TO RESTFRAME 


In order to compare the FP of the redshifted lens galaxies 
directly to the FP of the Coma cluster at 2 = 0.023, we 
calculate the effective surface brightness of the lens galaxies 
in restframe Johnson B band and restframe Gunn r band 
by interpolating between filters. 

In the following example we assume that the redshifted 
r band falls between the WFPC2 F555W (=V) and F814W 
(=7) filters and the effective surface brightness is determined 
in the NICM0S F160W (=77) reference filter. We assume a 
linear relation between the AB magnitudes of the redshifted 
r filter and the V and 7 filter, so that for the restframe r 
magnitude of a galaxy at redshift 2 we can write 

r z = V-a(V-1)-a(cv-ci) + cv — Cr + 2.51og(l + 2). (1) 


The constants c are the conversion constants between the 
standard Vega magnitudes and AB magnitudes 


c = 2.5 log 


( f 0 °°T(u)f^du\ 
{ f 0 °°T(v)dv ) 


- 48.60, 


( 2 ) 


with /J 6 ® 11 the fluxdensity of Vega, and T(v) the filter trans¬ 
mission. The transmission curves of the HST filters (including 
the CCD response) were obtained from STScI 3 , and those of 
the B and r passbands were obtained from Bessel (1990) and 
Thuan & Gunn (1976) respectively. We used the CALCPH0T 
task of the STSDAS package in IRAF to calculate the conver¬ 
sion constants for the HST filters. The conversion constants 
for the B and r passband follow from Bessel (1990) and Frei 
& Gunn (1994) respectively. The last term in 0 includes 
the broadening of the r band with redshift and makes the 
magnitude behave as if it is a flux, rather than a fluxdensity. 

Equation 0 relates the V — r z color of a lens galaxy 
to its observed V — I color. To determine a, we model the 


underlying populations and are less sensitive to outliers than stan¬ 
dard estimators (e.g. Beers, Flynn & Gebhardt 1990). 3 http://www.stsci.edu/instruments/observatory/cdbs/cdbs.html 
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Figure 1. Transformation to 
restframe r band by interpola¬ 
tion between the F555W (=V) 
and F814W (=1) filter. For 
four spectral types the V — I 
and V — r z colors are mod¬ 
eled (open circles). With the 
observed V — I color of the 
lens galaxy (filled circle) we in¬ 
terpolate between the modeled 
colors of the two nearest spec¬ 
tral types, to Slid the V — r z 
color of the lens galaxy. 


V—I and V—r z colors for four different spectral types, E/SO, 
Sbc, Scd and Im, using the spectral energy distribution of 
Coleman, Wu & Weedman (1980). This gives four estimates 
of a. Using the observed V — I color of the lens galaxy we 
interpolate linearly between the modeled colors of the two 
nearest spectral types (Fig. [3J. In this way we obtain the 
best-fit V — r z color and corresponding value of a for the 
lens galaxy. To estimate the uncertainty in the conversion 
constants we compare modeled colors with those predicted 
by Frei & Gunn (1994). We find that the differences are 
small and estimate the uncertainty at 0.02 magnitudes. 

The observed V — H color is then used to relate the 
effective surface brightness in the reference H filter to that 
of the restframe r band 

Me.r.-At e,H + (V-H)-(V-r z ). (3) 

Similar transformations are derived for each lensing system. 
In cases where observations are available in more than two 
passbands, we calculate H B ,r z for all filter combinations and 
make a selection based on the following criteria: the error 
in Ue,r z , the wavelength difference between the redshifted 
r band and the observed filter, and a preference for filter 
pairs enclosing the redshifted r band. The latter implies an 
interpolation between two filters, whereas otherwise we have 
to extrapolate. As byproducts of our procedure we find the 
SEDs that provide the best fits to the observed colors of the 
lens galaxies. For the restframe B and r band this yields 
respectively 20 (77%) and 19 (73%) lens galaxies that are 
best fitted by the E/SO spectral type, whereas the colors of 
the remaining lens galaxies are closest to the Sbc spectral 
type. 

Since r z in Q behaves like a flux, the effective surface 
brightness decreases as (l + z ) 4 with increasing redshift 4 . We 
correct the surface brightnesses for this cosmological dim¬ 
ming. For the 26 lens galaxies, the resulting values for the 
restframe Johnson B and Gunn r effective surface brightness 
He,B z and p e ,r z are given in Tabled together with the other 
two FP parameters; the central stellar velocity dispersion a c 
and effective radius r e . We have also included the redshift 
z of the lens galaxies and the velocity dispersions measured 
from stellar kinematics cr c * when available. 



Figure 2. Edge-on view of the FP in the restframe Gunn r band. 
The Coma cluster galaxies (small dots) and corresponding linear 
fit represent the local FP (j0rgensen, Franx and Kjaergaard 1995a, 
1995b, 1996). Lens galaxies for which the redshifts are known 
spectroscopically are indicated with circles. Filled circles corre¬ 
spond to lens galaxies for which the transformation to restframe 
(Section 0 was by interpolation between a pair of filters, and open 
circles if a modest extrapolation was needed. The open squares 
represent the 5 lens galaxies with a redshift that is estimated pho¬ 
tometrically (TableQJ. The transformation to restframe for these 
5 galaxies was by interpolation. For 7 lens galaxies, velocity dis¬ 
persions have been measured from stellar kinematics (see text for 
details). Asterisks show the position of these galaxies in the FP if 
the stellar velocity dispersion is used. The lens galaxies are offset 
from the Coma FP, as expected from evolution of their stellar 
populations. 

4 FP AND M/L EVOLUTION 

The Fundamental Plane has the form 

log r e = a log a c + f3 log 7 e + 7 (4) 

(Dressier et al. 1987; Djorgovski & Davis 1987), with r e in 
kpc, a c in kms -1 and p, e (= —2.51og/ e ) in magarcsec -2 . 
We adopt for the coefficients a and (3 the values derived 
by Jprgensen et al. (1996) for a sample of 225 early-type 
galaxies in nearby clusters. They found for the Johnson B 
band a = 1.20 ± 0.06 and (3 = —0.83 ± 0.02, and for the 
Gunn r band a = 1.24 ± 0.07 and (3 = -0.82 ± 0.02. 

We use the tabulated photometric and spectroscopic 
data of Jprgensen et al. (1995a, 1995b) to construct the FP 
of Coma. The edge-on projection of the Coma FP in the 
r band is shown in Fig. 01 (small dots). A linear fit to the 
coma FP yields an intercept of 9.50 ± 0.02 and 9.12 ± 0.02 
for the B and r band respectively. Large symbols show the 
lens galaxies. The lens galaxies show a large scatter, and are 
offset with respect to the FP of Coma. This relative differ¬ 
ence can be attributed to the evolution of the M/L ratios 
of galaxies, and the large scatter may in part be caused by 
the large range of redshifts spanned by the lens sample. 

4 The restframe r z magnitude 0 behaves therefore as a K- 
corrected magnitude. 
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As usual in FP evolution studies, we assume that early- 
type galaxies are a homologous family, i.e. that they are 
structurally similar. The total mass of a galaxy (including 
possible dark matter) is then proportional to an effective 
mass oc 1 j 2 r e . With the total luminosity proportional to I e rl, 
the FP relation implies that the effective mass-to-light ratio 
M/L oc M 02i r~ 002 in the r band (e.g. Treu et al. 2001). 
The tightness of the FP relation implies a low scatter in the 
M/L ratios of early-type galaxies of 23 % (Faber et al. 1987; 
Jprgensen, Franx & Kjaergaard 1996). Hence, the evolution 
of M/L can be well studied via the evolution of the FP. 

We assume that all early-type galaxies evolve in the 
same way, i.e. the coefficients a and /3 are independent of 
redshift and the same for cluster and held early-type galax¬ 
ies. Until now there is no convincing (observational) evi¬ 
dence against these assumptions, but this might change if 
more and deeper data become available (e.g. van Dokkum 
& Franx 1996; van Dokkum et al. 2001; Treu et al. 2002). If 
furthermore the effective radius and velocity dispersion do 
not change with redshift, the difference in FP intercept is 
due to a difference in surface brightness caused by luminos¬ 
ity evolution. As a result, the evolution of the intercept of 
the FP is proportional to the evolution of the mean M/L 
ratio. For each lens galaxy, the difference in FP intercept 
with respect to the Coma FP (Fig. |5J can be expressed as 
an offset in its M/L ratio relative to that of a galaxy of the 
same mass in Coma. Fig. [3 (two left panels) shows the de¬ 
pendence of the M/L offset on redshift, in the restframe r 
band. The M/L ratio of lens galaxies clearly shows a trend 
with redshift, with the highest redshift objects having the 
lowest M/L ratios. 

From a linear fit to the M/L ratios of the 26 lens galax¬ 
ies, we find an evolution rate dlog(M/L)/d« of —0.62 ±0.13 
and —0.47 ± 0.11 in restframe B and r band respectively. 
The intercept of the fit, A log M/L at 2 = 0, is —0.01 ± 0.08 
and —0.01 ±0.06 respectively, so that locally there is no sig¬ 
nificant offset in the mean M/L ratio of the lens and cluster 
galaxies. If this offset is forced to be zero, we find best-fit 
evolution rates for the B and r band of —0.63 ± 0.06 and 
—0.49 ±0.05, which are only slightly higher than for the un¬ 
constrained fit. (The lcr errors from the fit are smaller than 
for the unconstrained fit since only the slope is a fitting 
parameter, while the intercept is fixed to zero.) Restricting 
the fit to only the 21 lens galaxies with spectroscopic red- 
shifts, has no significant effect on the results. Similarly, we 
find no signihcant differences if we exclude the lens galaxies 
from the fit for which a modest extrapolation was needed 
in the transformation to restframe band; or if for the 7 lens 
galaxies for which the velocity dispersion <r c * measured from 
stellar kinematics is available, this value is used instead of 
the modeled a c from the lensing geometry. 

For cluster galaxies van Dokkum et al. (1998) find in 
the B band an evolution rate d \og(M/L)/dz = —0.49 ± 
0.05 This means on average a faster evolution for the lens 
galaxies, but the difference is only significant if the M/L 
ratio of lens galaxies is forced to coincide with that of cluster 
galaxies at z = 0. The value for the cluster evolution rate 
by van Dokkum et al. (1998) does not take into account 
possible effects due to morphological evolution. If early-type 
galaxies were transformed from late-type galaxies at modest 
redshifts (e.g. Dressier et al. 1997) the early-type galaxies 
that were already present at high redshift are only a subset of 


all progenitors of low redshift early-type galaxies. As a result 
of this “progenitor bias” (vDFOl) the formation redshift of 
morphologically selected cluster early-type galaxies may be 
overestimated. The lens galaxy sample is probably much less 
affected by progenitor bias, since they are selected on mass 
and not on morphology, but due to the merging of galaxies, 
a fraction of the progenitors may still not be accounted for. 
After applying the maximum correction for progenitor bias 
allowed by the data, vDFOl find a cluster evolution rate of 
dlog(M/ L)/dz = —0.56 ±0.05 in the B band, similar to the 
evolution rate for the lens galaxies. These results suggest 
that, if there is no progenitor bias, field galaxies may be 
younger than cluster galaxies, but that such an age difference 
becomes less significant if we correct for possible progenitor 
bias. 

Our results are consistent with determinations for field 
galaxies based on direct spectroscopic measurements of the 
velocity dispersions. For the B band, van Dokkum et al. 
(2001) arrive at an evolution rate d\og(M/L)/dz of —0.59± 
0.15, and Treu et al. (2002) find a value of —0.72lQ jg. Com¬ 
paring our results for the B band with those of R03, we find 
that the evolution rate obtained by R03 is on average slower, 
but the (lcr) confidence limits still overlap. For intercepts 
that are allowed to vary and forced to be zero, R03 find an 
evolution rate of —0.54±0.09 and —0.56±0.04 respectively. 

The error analysis and the transformation to restframe 
of R03 differs from ours (see also Section 0 , but we cannot 
further investigate the effects of these differences, as R03 do 
not give their resulting fundamental plane parameters. How¬ 
ever, if we add (in quadrature) an constant additional error 
to the observed uncertainties in the M/L ratios of the lens 
galaxies, the slope of the fit does decrease. Around the linear 
fit we measure a (biweight) scatter of 0.17 for the B band 
and 0.15 for the r band 5 , which in both cases is higher than 
the expected scatter from the observational errors of 0.12. If 
we now take for the constant additional scatter the difference 
(in quadrature) between the measured and expected scatter, 
we find for the unconstrained fit in the B band an evolution 
rate of —0.56 ±0.12, almost identical to the result of R03. If 
the offset is forced to zero the slope changes less, giving an 
evolution rate of —0.61 ± 0.05. Note that R03 rescale their 
input errors so that the best-fit model has a reduced of 
one. They assume that the additional scatter is due to un¬ 
derestimated errors in the data set, whereas we propose that 
internal population differences cause the additional scatter. 
Hence, we multiply our results with the square root of the 
reduced \ 2 to reflect this aspect, but we do not change the 
uncertainties on the input data. 

The additional scatter implies that the M/L ratios of 
the lens galaxies are not well fitted with a single evolution 
rate, whereas for the cluster galaxies the fit is very good (e.g. 
van Dokkum et al. 1998). This may be due to a significant 
spread in stellar population ages among field galaxies, which 
induces an additional scatter measured in the evolution rate 
dlog(M/ L )/dz of the lens galaxies In the next Section, we 
relate the M/L evolution to stellar population ages by fit¬ 
ting simple single burst models. In Section|S] we then study 


5 The scatter in the B band is probably higher than in the r 
band due to the fact that the B band is more sensitive to recent 
star formation. 


the restframe colors of lens galaxies and their evolution. To 
find out whether there is a significant age spread, we also 
investigate scatter in color and if it is correlated with the 
scatter in M/L ratio. 


5 STELLAR POPULATION AGES 

The evolution of the M/L ratio depends on the age of the 
stellar population. A stellar population which formed at 
low redshift will evolve faster than a population formed at 
high redshift: the luminosity of a young population becomes 
rapidly fainter when short-lived massive and bright stars dis¬ 
appear, whereas the dimming of the light is more gradual 
for an old population dominated by low mass stars. We esti¬ 
mate the stellar population ages of lens galaxies and the age 
difference between lens and cluster galaxies by fitting simple 
single burst models to the M/L evolution. 

The M/L ratio of a single burst stellar population with 
fixed mass evolves as 

M/L oc (t — t*) K , (5) 

with t* the stellar formation time, corresponding to a stellar 
formation redshift z* (e.g. Tinsley & Gunn 1976). The co¬ 
efficient k depends on the Initial Mass Function (IMF) and 
metallicity, and also on the passband in which the lumi¬ 
nosity is measured. For a normal IMF with Salpeter (1955) 
slope and Solar metallicity, kb ~ 0.93 and K r ~ 0.78 for 
the restframe B and r bands (Worthey 1994). Note that the 
predicted evolution is independent of Ho because the age 
dependence of the M/L ratio is a power-law. 

We first investigate if we can fit the observed M/L evo¬ 
lution of the lens galaxies by a single burst evolution model 
with the same stellar formation redshift z* as has been de¬ 
rived previously for cluster galaxies (the null-hypothesis). In 
Section 5.2 and 5.3 we investigate the range of z* allowed 
by the data. 

5.1 Can cluster galaxies and lens galaxies have 
the same age? 

Before investigating more complex models we consider the 
case that lens galaxies and cluster galaxies have the same lu¬ 
minosity weighted stellar age. To determine the mean stellar 
formation redshift of cluster early-type galaxies (zj*) we fit 
a single burst model 0 to previously published restframe B 
band data for the clusters CL1358+62 (Kelson et al. 1997), 
CL0024+16 (van Dokkum & Franx 1996), MS2053+03 (Kel¬ 
son et al. 1997) and MS1054-03 (van Dokkum et al. 1998) at 
redshifts 0.33, 0.39, 0.58 and 0.83, respectively, normalized 
with respect to the average M/Lb ratio of the local Coma 
cluster at z = 0.023 (Jprgensen et al. 1996). 

To perform such a single burst fit, we minimize 

2 ^ ( M/L mo d,i - M/Lobs,; \ 2 

X -g(- „(M/W) -) • (6) 

with n the number of M/L observations used as constraints, 
M/L a bs,i the *th observation, M/L mo d,i the corresponding 
single burst prediction and a(M/L 0 h s ,i) the uncertainty or 
error in this observation. To determine the confidence levels, 
we calculate the difference in X' 2 between a model and the 
overall minimum, A\ 2 = \ 2 ~ Xmim to which the usual 
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Gaussian confidence levels can be assigned (e.g. Press et al. 
1999). 

The resulting best-fit mean formation redshift (zj f ) = 
4-3 ^i‘ 2- This direct fit does not take into account possible 
progenitor bias. After applying the maximum correction for 
progenitor bias allowed by the data, vDFOl find (zj ! ) = 
2-0^0 2 f° r duster galaxies. Therefore, we consider 2.0 < 
zf < 4.3 as a plausible range for the mean star formation 
epoch of cluster early-type galaxies. 

We fitted single burst models with this range of for¬ 
mation redshifts to the restframe B and r band data of 
the lens galaxies. The fits improve towards higher forma¬ 
tion epoch, but even those with maximum progenitor bias 
correction are rejected with nearly 100 %. We can allow for 
an offset between M/L ratios of lens galaxies and those of 
cluster galaxies. We will describe this offset with the value 
of A log M/L at z = 0, which we denote by b. For a stellar 
formation epoch ranging from 4.3 to 2.0 (no to maximum 
progenitor bias correction) the offset b varies from —0.12 to 
—0.06 for the B band, and from —0.08 to —0.03 for the r 
band. These models with a systematic offset fit the lens data 
better, but are still rejected with > 99.9 % confidence. 

5.2 Model A: unconstrained fit 

We assume that the M/L ratios of field galaxies are inde¬ 
pendent of those of cluster galaxies, i.e. a cluster galaxy of 
a given age can have a very different M/L ratio than a field 
galaxy of the same age. This may be the case if, e.g. the 
metallicities of field and cluster galaxies are different at a 
given mass. In addition to the stellar formation redshift z*, 
we also have the normalization of the single burst model as 
free parameter. We describe this second parameter with b 
(see Section 15.11 1. 

For a range of z* and b values, we fit single burst mod¬ 
els the observed M/L values of the redshifted lens galaxies. 
For the r band, the 1,2 and 3 cr limits on z* and b are 
shown the in the upper-right panel of Fig. 0 with the min¬ 
imum indicated by a cross. The best-fit values of z* and b 
are l.SjfJ'g and —0.03 ± 0.09 for the B band, and 1.9ftJ’e 
and —0.01 ± 0.07 for the r band (Tabic |2J. Note that the 
given uncertainties are lcr errors unless noted. The M/L r 
evolution that corresponds to the best-fit values is shown in 
the upper-left panel of FigHTl (dashed line). Since for both B 
and r band the model parameter b is not significantly dif¬ 
ferent from zero, we obtain similar best-fit values for z* if 
we normalize the single burst model such that 6 = 0, i.e. if 
we assume that locally the average M/L ratio of field and 
cluster galaxies is the same. 

5.3 Model B: simultaneous fit to lens and cluster 
galaxies 

Here we assume that the stellar populations of field and 
cluster galaxies evolve in the same way. The stellar popu¬ 
lations of the lens galaxies may form at a different redshift 
than those of the cluster galaxies, but galaxies of a given age 
have identical M/L ratios. For the formation redshift of the 
stars in cluster galaxies we use the values that we obtained 
in Section rm z? = 4.3 if not corrected for progenitor bias 
and z? = 2.0 after correction for maximum progenitor bias. 
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Figure 3. Upper-Left: The evolution of the M/L offset in the restframe Gunn r band. The small dot is the average M/L of the Coma 
galaxies. The lens galaxies are indicated with the same symbols as in Fig.\2\ The solid line represents the single burst evolution of a stellar 
population formed at the mean cluster formation redshift of z+ = 4.3. If we correct for progenitor bias, the latter formation redshift will 
be lower, and hence the evolution steeper. This is indicated by the shaded region, with z* = 2.0 for maximum progenitor bias correction 
(vdFOl). The dashed line is the best fit model for the lens galaxies, with z+ = 1.9 and zeropoint b = A log M/L rz (z = 0) = —0.01. 
Upper-Right: Confidence levels (1, 2 and 3a) for combinations of z* and the zeropoint b when they are both allowed to vary. Lower-Left: 
The solid line shows the evolution for a stellar population formed at z* = 2.0, which is the mean cluster value after maximum progenitor 
bias correction. The dashed line is the model for a stellar population formed at the best-fitting star formation redshift of z+ = 1.8, if the 
lens and cluster galaxies are fitted in a self-consistent way (model B). To illustrate the effect of changing z* we also show a model with 
z* = 1 (dotted line). Lower-Right: Because cluster and lens galaxies are modeled self-consistently, the stellar formation redshift of the 
lens galaxies depends on that of the clusters z+ l . The best-fit values vary from 1.8 to 2.3 (indicated by arrows) if z% 1 is increased from 2.0 
to 4.3, the values for maximum and no progenitor bias correction. The errors on the best-fit parameters from the confidence levels are 
multiplied with the square root of the reduced \ 2 of the fit, to reflect the measured additional scatter compared to the expected scatter 
from the observational errors (see text and Tabled for resulting uncertainties on b and z*). 
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model 


B z band 

r z band 

A 


1.8 (-0.5/+1.4) 

1.9 (-0.6/+1.9) 

B (with z/ 1 

= 2.0) 

1.7 (-0.1/+0.2) 

1.8 (-0.2/+0.2) 

B (with z// 

= 4.3) 

2.1 (-0.2/+0.3) 

2.3 (-0.3/+0.3) 


Table 2. The mean formation redshift (z*) of a single burst stellar 
population for lens galaxies (la errors). Models as in Fig. 0 (see 
text for further details). 


In both cases we determine the constraints on z* for the lens 
galaxies. Note that z* and b are coupled in this model, be¬ 
cause of the constraint that lens galaxies with the same age 
as cluster galaxies have identical M/L ratios. For the r band 
the resulting constraints on z* (and hence b) are shown in 
the lower-right panel of Fig. [3 with the best fit values in¬ 
dicated by arrows. From maximum to no progenitor bias 
correction, we find z* ranging from 1.6 to 2.4 (la) for the 
B band, and from 1.6 to 2.6 (la) for the r band (Tabled- 
For the case of maximum progenitor bias the M/L evolution 
that follows from the best-fit value z* = 1.8 in the r band, 
is shown in the lower-left panel of Fig. 0 (dashed line). 

5.4 Summary of results 

We have demonstrated that the M/L evolution of lens galax¬ 
ies cannot be fitted with models that provide good fits to 
cluster galaxies. The fit clearly improves if lens galaxies are 
allowed to be systematically offset from cluster galaxies due 
to metallicity differences or other systematic effects. 

If we consider both this offset and the stellar formation 
epoch of the lens galaxies as free parameters, the best-fit 
single burst model in the B and r band are consistent. The 
resulting offset is not significantly different from zero, and 
although the best-fit stellar formation redshift implies for 
lens galaxies on average a younger stellar population than for 
cluster galaxies, the resulting (la) range of 1.3 < (z*) < 3.8 
is not significantly different from the formation epoch for 
cluster galaxies of 1.8 < (z//) < 8.0 for the full range from 
no to maximum progenitor bias correction. 

If we impose the constraint that galaxies of the same age 
have the same M/L ratio irrespective of their environment, 
we find stellar formation redshifts for the lens galaxies of 

1.6 < (z*) < 2.0 if the stars in cluster galaxies formed at 
z / l = 2.0 (maximum progenitor bias), and 1.9 < (z*) < 

2.6 if the stars in cluster galaxies formed at z= 4.3 (no 
progenitor bias). In the latter case the stellar populations 
of the lens galaxies are significant younger (10-15% at the 
present epoch) than those of the cluster galaxies. 

If the local M/L offset between lens and cluster galax¬ 
ies is allowed to vary, R03 find for the stellar formation red¬ 
shift of the lens galaxies a (la) range of 2.0 < (z*) < 3.6 
from their B band analysis. Although we find a range of 
1.3 < (z*) < 3.2 (model A, B band) which implies on aver¬ 
age a somewhat younger stellar population, the results are 
consistent. R03 conclude that the M/L evolution rates they 
measure favor old stellar populations for the lens galaxies 
with a mean formation redshift (z*) > 1.8 at a 2a confidence 
level. We find a lower 2a confidence limit of (z*) > 1.2. 

To test the hypothesis whether the M/L evolution of 
the lens galaxies can be fitted with a single burst model, 


we calculate the reduced y 2 . We also compare the measured 
scatter around the fit with the expected scatter from the 
uncertainties in the M/L ratios of the lens galaxies. If we al¬ 
low more freedom in the single burst models, the fit improves 
(lower reduced \ 2 ) and the measured scatter decreases. How¬ 
ever, even the over-all best-fit single burst model is rejected 
with > 99 % confidence and the measured scatter of 0.17 and 
0.15 in B and r band is significant higher than the expected 
scatter of 0.12 from the observational errors. To establish 
whether the additional scatter is due to a significant spread 
in ages among field galaxies, we study in the next Section 
the (restframe) colors of the lens galaxies, and investigate 
if there is also a significant scatter in the color and if it is 
correlated with the scatter in their M/L ratios. 


6 COLORS 

If the stellar populations of lens galaxies are on average 
younger than those of cluster galaxies (Section^ we expect 
their colors to evolve more rapidly and to be on average 
bluer than those of cluster galaxies. For single burst stellar 
populations the B — r color evolves as 

B — r = 2.5(kb — Kr) log(f — f*) + c, (7) 

with k,b — Kr ~ 0.15 (Worthey 1994) and c a normalization 
constant. 

The left panel of Fig.0]shows the restframe B — r colors 
of the lens galaxies versus redshift. For the lens galaxies 
with open circles, a modest extrapolation was required in 
the transformation to either the restframe B or r band. The 
red outlier (at z ~ 0.4) is the lens BRI0952-0115, for which 
the observed R — H (F675W-F160W) color is significantly 
redder than the modeled E/S0 color, and for which we had 
to extrapolate in the transformation to restframe B band. 
The single burst evolution of a stellar population formed at 
the mean cluster formation epoch of zf = 2.0 (maximum 
progenitor bias), is drawn (solid line) through the averaged 
color of the coma cluster galaxies (small dot). A decline in 
the colors of the lens galaxies with increasing redshift is 
visible, in spite of the large scatter. If we fit single burst 
models as in Section 0 (model A and B), we find on average 
a younger stellar formation epoch for the lens galaxies, but 
the difference with the cluster galaxies is never significant. In 
case the stellar populations of field and cluster galaxies are 
assumed to evolve in the same way (model B), the evolution 
model with the best-fit stellar formation redshift of z* = 1.7 
is shown in the left panel of Fig. 0 (dashed line). Similarly as 
for the single burst model fits to the M/L evolution, the fits 
to the color evolution are rejected with > 99 % confidence, 
and the measured scatter is in all cases significantly higher 
than the expected scatter from the observational errors in 
the restframe colors of the lens galaxies. 

To further investigate the scatter in color, we subtract 
the predicted colors of cluster galaxies with zf = 2.0, from 
the restframe colors of the lens galaxies. Note that this 
choice of formation redshift corresponds to the minimum 
age difference between field and cluster galaxies. In the right 
panel of Fig. 0 we show the resulting residual colors plotted 
against velocity dispersion, with a linear fit to the Coma 
cluster galaxies. The lens galaxies show a large spread in 
their reduced colors, and are on average bluer by ~ 0.1 mag 
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Figure 4. Left: Evolution of the restframe B — r color. The symbols correspond to those in Fig. ^ with the colors of the coma cluster 
galaxies averaged (small dot). The solid line shows the evolution for a stellar population formed at 2 * = 2.0, which is the mean cluster 
value after maximum progenitor bias correction. The dashed line is the model for a stellar population formed at the best-fitting star 
formation redshift of z* = 1.7, if the lens and cluster galaxies are fitted in a self-consistent way. To illustrate the effect of changing 2 * we 
also show a model with 2 * = 1 (dotted line). For redshifts beyond z ~ 0.5 the scatter in color increases and becomes significant higher 
than expected from the observational errors. Right: The restframe B — r colors of lens galaxies, after subtracting the fiducial model for 
cluster galaxies with z ( / = 2.0, plotted versus velocity dispersion. The solid line is a linear fit to the Coma cluster galaxies. On average 
the colors of the lens galaxies are bluer than those of the cluster galaxies, consistent with a younger average stellar population. The 
scatter in the reduced colors is larger than expected from the observational errors. 


than the cluster galaxies. The bluer average color is quali¬ 
tatively consistent with the on average younger ages of lens 
galaxies derived from the single burst model fits above and 
from the fits to the M/L evolution (Section 15.41 . For the 
total sample of 26 lens galaxies, we measure a scatter in 
the residual colors of 0.22, that is significant higher than 
the expected scatter of 0.18. From the left panel of Fig. 2] 
we observe that the color scatter increases beyond z ~ 0.5. 
For the 15 lens galaxies with a > 0.5 we measure a scat¬ 
ter of 0.24, whereas for the lens galaxies with z < 0.5 the 
measured scatter is only 0.16 and nearly identical to the ex¬ 
pected scatter of 0.15 from the errors in the colors of these 
galaxies. 

The additional measured scatter in the residual colors 
of the lens galaxies may indicate a significant spread in the 
ages of the stellar populations of the lens galaxies. We test 
whether the color spread is caused by a spread in ages or 
other effects, by investigating whether the residual colors 
correlate with the residual M/L ratios. In Fig. 0 we plot 
the color of the lens galaxies against the logarithm of their 
M/L ratio. For both quantities the expected evolution of 
a stellar population formed at z// = 2.0 was subtracted; 
hence a galaxy with the average color and M/L ratio of 
galaxies in Coma would be located at the origin. From Eq. 
0 and 0 it follows that a linear relation is expected in Fig. 
[S] with a (time-independent) slope of 2.5 (kb — K r )/K>B ~ 
0.40 (solid line). The data appear to be correlated in the 
expected sense, albeit with large scatter. To test whether the 
correlation is significant we used the Spearman’s rank-order 
correlation coefficient rs- We find that rs = 0.47, so that 
with N = 26 lenses, we can reject the hypothesis that the 



Figure 5. The B — r color versus M/L r ratio, after subtracting 
the fiducial cluster galaxy evolution model with z+ l = 2.0. The 
symbols are the same as in Fig. [3 with the average of Coma 
situated in the origin. The dashed line is a linear fit to the data; 
the solid line shows the expected correlation between color and 
M/L ratio due to age variations, with slope 2.5 (kb ~ Kt)/kb ~ 
0.40. Although the slope is uncertain, the correlation between 
color and M/L ratio is significant at the 95 % level. Hence we 
interpret the intrinsic scatter in the colors of the lens galaxies as 
a stellar population effect, probably caused by a spread in their 
ages. 
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two quantities are uncorrelated with > 95 % confidence. The 
only viable explanation for the correlation is age variation. 
If the correlation would be caused by metallicity variations, 
(some) field galaxies would be much less metal rich than 
cluster galaxies, opposite to the result by Kuntschner et al. 
( 2002 ). 

The colors and M/L ratios of the bluest lens galaxies 
are best fitted with stellar formation redshifts as low as 2 * ~ 
1. About half the lens galaxies are consistent with an old 
cluster-like stellar population with stellar formation redshift 
Zi, > 2. If galaxies form in a sequence of bursts, formation 
redshifts are indicative of the last prominent epoch of star 
formation. The galaxies can have older underlying stellar 
populations. 


7 SUMMARY AND CONCLUSIONS 

We studied the evolution of the M/L ratios of lens galaxies 
in the restframe Johnson B and Gunn r bands. For an flat 
cosmology with Qm = 0.3 and LIa = 0.7, we obtained an 
evolution rate dlog(M/ L)/dz of —0.62 ± 0.13 in restframe 
B and —0.47 ±0.11 in restframe r. Due to differences in the 
determination of the FP parameters and the corresponding 
errors, the evolution rate obtained by R03 of —0.54 ±0.09 is 
slightly slower but not significant different. Our results are 
consistent with determinations for field galaxies based on di¬ 
rect spectroscopic measurements of the velocity dispersions. 
For the B band, van Dokkum et al. (2001) arrive at an evo¬ 
lution rate of —0.59 ± 0.15, and Treu et al. (2002) find a 
value of —0.72^Q')g. The weighted mean of these results and 
our result yields an evolution rate for field early-type galax¬ 
ies of —0.64 ± 0.06 in the B band. For cluster galaxies van 
Dokkum et al. (1998) and vDFOl find an evolution rate be¬ 
tween dlog(M/Ls)/d« = —0.49±0.05 and —0.56±0.05, for 
minimum and maximum progenitor bias respectively. The 
latter value is not significantly different from the average 
evolution rate for the field galaxies. 

We further investigated the M/L evolution in Section 
0 where we related it to stellar population ages by fitting 
simple single burst models for a Salpeter (1955) IMF and 
Solar metallicity. The M/L evolution of cluster galaxies is 
well approximated by a mean single burst formation redshift 
of z/ 1 = 4.3j^'2, and zf = 2.01q 2 a ft er maximum progenitor 
bias correction (vDFOl). 

We first tested if these cluster models could fit the M/L 
evolution of the lens galaxies. This is not the case, but the 
fits improve if there is a systematic offset between lens and 
cluster galaxies of ~ —0.1 in A log M/L Such an offset could 
be caused by, e.g. metallicity variations, systematic differ¬ 
ences in the velocity dispersions due to the different mea¬ 
surement techniques, or other effects. It is interesting to note 
that hierarchical models have predicted such a constant off¬ 
set with redshift between cluster and field galaxies (see van 
Dokkum et al. 2001). However, similar to van Dokkum et al. 
(2001), we conclude that the observed offset is much smaller 
than the predicted offset of A log M/Lb ~ —0.26. 

If we next allow the stellar formation redshift 2* also to 
vary, we find for the best-fit single burst model that the off¬ 
set is not significant from zero, together with a mean stellar 
formation epoch ( 2 *) for the lens galaxies of 1.8lJ'g in the 
B band and 1.9jtjg in the r band. This means that on aver¬ 


age the stellar populations of the lens galaxies are younger, 
but the difference with the cluster galaxies is not signifi¬ 
cant. However, if we impose the constraint that galaxies of 
the same age have the same M/L ratio irrespective of their 
environment, we find without correction for progenitor bias 
that the stellar populations of the lens galaxies are signifi¬ 
cant younger (10-15% at the present epoch) than those of 
the cluster galaxies. In the case of maximum progenitor bias 
the average difference (~ 5%) is not anymore significant. 

From their analysis R03 obtain a (2a) lower limit for the 
mean stellar formation epoch of lens galaxies of ( 2 *) > 1.8, 
whereas we find that lower stellar formation redshifts are 
allowed, with ( 2 *) > 1.2 as a 2a lower limit. Nevertheless, 
these results disagree with the prediction of semi-analytical 
hierarchical galaxy formation models (e.g. Kauffmann 1996; 
Kauffmann & Chariot 1998; Diaferio et al. 2001) that early- 
type field galaxies in general have very late star formation 
with 2 * < 1. 

Whereas the M/L evolution of cluster galaxies is well 
approximated by a single burst evolution model, we found 
that this is not the case for the lens galaxies. All single burst 
models are rejected with >99% based on the reduced x 2 , 
and the measured scatter is in all cases significant higher 
than the expected scatter from the observational errors. The 
additional scatter is most likely caused by differences in the 
stellar population of the lens galaxies. 

To investigate whether there is a significant spread in 
ages, we studied in Section|3the restframe colors of the lens 
galaxies and their evolution. The colors of the lens galax¬ 
ies are on average bluer than those of the cluster galaxies, 
consistent with on a younger average stellar population. We 
found that the measured scatter in the single burst fits to 
the colors is much larger than the expected scatter from the 
observational errors. Moreover, we showed that there is a 
significant correlation between the colors and the M/Lb ra¬ 
tios of the lens galaxies. We interpret this as evidence for a 
significant spread in the stellar population ages of the lens 
galaxies. 

Whereas about half of the lens galaxies are consistent 
with old cluster-like stellar populations, the bluest galaxies 
are best fit by single burst models with young stellar forma¬ 
tion redshifts 2 * ~ 1. For the seven blue lens galaxies with 
residual colors below —0.2 mag in the right panel of Fig. 
Q] we found (as a byproduct of our transformation to rest¬ 
frame B band) that only two of them are best fitted by the 
E/S0 spectral type, i.e. ~ 28%, whereas for the total sam¬ 
ple we found 77 %. For the r band even six of the seven blue 
galaxies are best fitted by the Scd spectral type, instead of 
the ES/0 spectral type. Three of the seven blue galaxies in¬ 
deed show (some) deviations from early-type morphology. 
FBQ0951±2635 is an edge-on disk galaxy, SBS1520±530 
is slightly irregular and B1608±656 is an apparently dusty 
galaxy with star forming regions. However, a homogeneous 
sample of quasar subtracted (NICM0S) images of the lens 
galaxies is needed to do a more detailed and systematic 
study of their morphologies. Such a study will be valuable 
to determine the cause of apparently younger populations 
in a fraction of the lens galaxies. 
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